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Abstract 



Charge Coupled Devices (CCDs) have been successfully used in several high energy physics 
experiments over the past two decades. Their high spatial resolution and thin sensitive layers 
make them an excellent tool for studying short-lived particles. The Linear Collider Flavour 
^ Identification (LCFI) collaboration is developing Column-Parallel CCDs (CPCCDs) for the 
vertex detector of a future Linear Collider. The CPCCDs can be read out many times faster 
than standard CCDs, significantly increasing their operating speed. An Analytic Model has 
been developed for the determination of the charge transfer inefficiency (CTI) of a CPCCD. 
The CTI values determined with the Analytic Model agree largely with those from a full 
TCAD simulation. The Analytic Model allows efficient study of the variation of the CTI on 
parameters like readout frequency, operating temperature and occupancy. 
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Abstract — Charge Coupled Devices (CCDs) have been suc- 
cessfully used in several high energy physics experiments over 
the past two decades. Their high spatial resolution and thin 
sensitive layers make them an excellent tool for studying short- 
lived particles. The Linear Collider Flavour Identification (LCFI) 
collaboration is developing Column-Parallel CCDs (CPCCDs) for 
the vertex detector of a future Linear Collider. The CPCCDs can 
be read out many times faster than standard CCDs, significantly 
increasing their operating speed. An Analytic Model has been de- 
veloped for the determination of the charge transfer inefficiency 
(CTI) of a CPCCD. The CTI values determined with the Analytic 
Model agree largely with those from a full TCAD simulation. The 
Analytic Model allows efficient study of the variation of the CTI 
on parameters like readout frequency, operating temperature and 
occupancy. 



I. Introduction 

Charge transfer inefficiency (CTI) is an important aspect in 
the CCD development for operation in High Energy Physics 
colliders [l]-[3]. The LCFI collaboration has been developing 
new CCD chips and testing them for about 10 years [l]-[5]. 
Recently the focus of the simulations has been on CCDs with 
column parallel readout (CPCCD). Full TCAD [6] simulations 
for a CPCCD were performed for different readout frequencies 
and operating temperatures [7]-[9]. An example of the CTI 
temperature dependence is shown in Fig. [T] (from [9]). Full 
TCAD simulations are very CPU intensive. This has already 
been noted for the CCD simulations with a sequential read- 
out [10]— [12]. The CTI depends on many parameters, such as 
readout frequency and operating temperature. Some parame- 
ters are related to the trap characteristics like trap energy level, 
capture cross-section and trap concentration (density). Other 
factors are also relevant, such as the occupancy of the pixels 
(hits). It is well known that analytic charge transfer models 
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Fig. 1. CTI determined with TCAD simulations as a function of temperature 
for a two-phase CPCCD for two traps, 0.17 eV and 0.44 eV, with a concen- 
tration of 10 12 cm -3 and 1% hit (pixel) occupancy at readout frequencies 
10, 25 and 50 MHz. 



can be used to study the CTI dependence on readout frequency 
and operating temperature [13]— [15]. For a comparison with 
full TCAD CTI simulation results, as shown in Fig. [T] we 
have developed Analytic Models for the CPCCD [7]-[9]. 
The further development of these Analytic Models leads also 
to better understanding of the relevant parameters in order 
to reduce the CTI in future CPCCD prototypes. This paper 
addresses the inclusion of signal shape and clock voltage 
amplitude which leads to an improved Analytic Model for a 
CPCCD. The CTI obtained from the improved Analytic Model 
is compared with results from full TCAD simulations. 



II. Analytic Model for CTI determination 

The Analytic Models [7]-[9] describe the different steps 
in the charge transfer process and the amount of the trapped 
charge with respect to the charge cloud in transfer. Figure [2] 
shows the consecutive charge transfer for a two-phase CPCCD 
in one pixel (2 nodes). Following the treatment by Kim [16], 
based on earlier work by Shockley, Read and Hall [17] a 
defect at an energy below the bottom of the conduction band 
is considered. Our model considers one single energy level 
and includes the emission time r e , and capture time r c , in the 



differential equation 



dr 
dt 



r 



(1) 



where r is the fraction of filled traps. Initially the fraction 
of filled traps is r(0). At stage A the signal charge packet 
arrives and interacts with traps under node 1 during time t\. 
This interaction leads to the capture and emission process. By 
resolving the differential equation (1), the fraction of filled 
traps t\a under node 1 during time t\ (when the signal packet 
is present) is given by 



r 1A (OHK0)--]exp(-^) + 



(2) 



where r s = r c r e /(r c + r e ). 

At stage B charge moves to the next node and interacts 
with traps during time t 2 under this node. During this time 
electrons emitted from node 1 join the signal charge packet in 
the second node. Thus, the fraction of filled traps t\b under 
node 1 during time t\ in the presence of the signal packet is 
given by 



riBfo) = ri A (h)ex-p(- — ). 



(3) 



At the same stage B, t 2 b is defined as the fraction of filled 
traps under node 2 during time t 2 , thus, 



r 2B (t 2 ) = [r(0) 
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When the signal charge moves to the first node of the next 
pixel, stage C, electrons emitted during time t\ can join the 
signal present at this node and the fraction of filled traps r 2 c 
under node 2 during time t\ is given by 



r2c(h) = r 2B (t 2 ) exp(- — ). 



(5) 



The CTI is defined by the ratio of the charge loss under each 
node to the signal charge density n s , thus, 



CTI 
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(6) 



where N t is the trap concentration, and r(0) = exp(— t w /r e ) 
which is determined by considering the fact that initially all 
traps are filled and electrons are emitted during the waiting 
time t w between two signal charge packets. For the case ti = 
t 2 = t, the combination of the previous equations leads to 
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Fig. 2. Diagram of consecutive transfer of charge in a two-phase CCD. The 
diagram shows the charge transfer at different stages in time and space. 
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Fig. 3. CTI values from an Analytic Model as a function of temperature in a 
two-phase CPCCD for the two traps, 0.17 eV and 0.44 eV with a concentration 
of 10 12 cm -3 and 1% hit (pixel) occupancy at readout frequencies 10, 25 
and 50 MHz. 



III. Analytic Model CTI results 

The CTI dependence on readout frequency and operating 
temperature has been explored using an Analytic Model based 
on Eq. (7). Figure [3] shows the CTI results from the Ana- 
lytic Model at different frequencies for temperatures between 
100 K and 550 K. The CTI increases as the readout frequency 
decreases. For higher readout frequencies there is less time 
to trap the passing signal, thus the CTI is reduced. At high 
temperatures the emission time is so short that the trapped 
charges can rejoin the passing signal. 

IV. Comparison between full TCAD simulation 
and Analytic Model regarding signal shape 

EFFECT 

The signal charge profile varies in the signal cloud as 
illustrated in the upper part of Fig. [4] The signal packet does 
not have well defined boundaries and the charge concentration 
decreases gradually from the centre of the signal packet. 
Therefore, the signal packet will interact with a varying 
fraction of the traps within the pixel and this affects the CTI 
determination. The implementation of a more realistic signal 
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Fig. 4. Diagram of the signal packet in a potential well and effect of the 
signal shape on the expected CTI. Upper part: the density of a signal packet in 
a potential well decreases gradually from the centre of the packet. Lower part: 
expected scenario of the CTI dependence on the shape of the signal packet 
(small CTI for a narrow shape as shown on the left-hand side and large CTI 
for a wider shape as shown on the right-hand side). 

shape into the Analytic Model is expected to improve the 
agreement with the full TCAD simulation. Figure [5] shows 
the profile of the signal charge under the node from a full 
TCAD simulation. Two-dimensional and one-dimensional sig- 
nal charge density profiles are extracted as shown in Figs. [6] 
and [7] respectively. 
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Fig. 5. Signal charge density in transit from a TCAD simulation of a CPCCD. 
The plot shows the charge packet located under a node at a depth of about 0.5 
microns. One pixel is located between x = 10 and 30 microns. The arrow 
indicates the direction of the transfer. 

Figures [8] and [9] show the CTI dependence on the signal 
charge profile for the 0.17 eV and 0.44 eV traps at 50 MHz. 
These figures also show the CTI values for different signal 
shapes in comparison with the full TCAD simulation. The 
CTI is reduced as the width of the potential well becomes 
smaller. This behaviour is expected, as illustrated in the lower 
part of Fig. [4] The CTI values calculated with the Analytic 
Model including the signal charge profile agree better with 
the full TCAD simulation results. The relatively shallow traps 
(0.17 eV) are more affected by the signal charge shape than 
the deeper ones (0.44 eV). The inclusion of the approximate 
signal shape in the Analytic Model reduces the CTI value in 
the peak region by about 10 to 20% compared to assuming a 
square- shape signal. 



Fig. 6. Two-dimensional signal charge density extracted from the charge 
packet under one node using a full TCAD simulation. 
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Fig. 7. One-dimensional signal charge density extracted from the charge 
packet under one node using a full TCAD simulation at a depth of 0.5 microns. 
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Fig. 8. CTI from Analytic Model (AM) including the shape of the signal 
packet as a function of temperature for 0.17 eV traps with a concentration of 
10 12 cm -3 and 1% hit (pixel) occupancy at 50 MHz readout frequency in 
comparison with full TCAD simulation results. Three different signal shapes 
are compared with the full TCAD simulation. The CTI calculation with the 
Analytic Model using the signal shape extracted from a full TCAD simulation 
agrees better with the full TCAD simulation results than those from the 
Analytic Model using a square-shape signal as assumed previously. 
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Fig. 9. CTI from Analytic Model (AM) including the shape of the signal 
packet as a function of temperature for 0.44 eV traps with a concentration of 
10 12 cm -3 and 1% hit (pixel) occupancy at 50 MHz readout frequency in 
comparison with full TCAD simulation results. Three different signal shapes 
are compared with the full TCAD simulation. For the 0.44 eV traps the 
inclusion of the signal shape in the Analytic Model has only a small effect 
to improve the agreement with the full TCAD simulation. 



Fig. 11. CTI from Analytic Model (AM) including clock voltage effects as a 
function of temperature for 0.17 eV traps with a concentration of 10 12 cm -3 
and 1% hit (pixel) occupancy at 50 MHz readout frequency in comparison 
with full TCAD simulation results. Two different clock voltages (%) are 
shown. 
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Fig. 10. Diagram of clock voltages in a two-phase CPCCD. V\ and V2 
(dashed lines) show the applied voltages under node 1 and 2, respectively. 
The solid line shows the difference between the two applied voltages. Vb is 
the barrier potential (horizontal dashed line). T\ and T3 are time periods with 
no charge transfer and T2 is the time period with charge transfer. 



V. Comparison between full TCAD simulation and 
Analytic Model regarding clock voltage effect 

In this study the effect of different clock voltage amplitudes 
on CTI values are investigated. A sine-form voltage is applied 
to consecutive nodes as shown in Fig. 
variables are defined: 



10 The following 



voltage applied to a first node of a pixel, 

voltage applied to a second node of a pixel, 

potential barrier created between two successive gates by 

the doping profile, 

time interval where |Vi — V2I < Vb, 

time interval where | V\ — V2 \ > Vb- 



The signal is not transferred until the absolute difference 
between the two clock voltages V\ and V2 reaches the potential 
barrier created between two consecutive nodes. This affects 
the CTI determination and it is now included in the Analytic 
Model. The time intervals Ti and T3 are defined by the 
intersection point between the | V\ — V2 1 curve and the barrier 
potential Vb (horizontal dashed line): Vm sm(wt) = Vb, thus, 
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x sin 
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where Vm is the amplitude of the clock voltage. The CTI 
determined with the Analytic Model including the clock 



voltage effect is shown in Figs. 11 and 12 for 0.17 eV and 
0.44 eV traps, respectively. These results are compared to full 
TCAD simulations. Two different clock voltage amplitudes are 
used to illustrate the effect of the clock voltage. It is noted that 
the CTI decrease occurs only for temperatures above the CTI 
peak position. 

In addition, the effect of the clock voltage amplitude on 
the CTI is studied for a large variation of amplitudes. The 
CTI decreases as the amplitude increases until it saturates and 
no further decrease can be observed. This result is shown in 



Fig. 13 for two examples, 0.17 eV traps at a temperature of 
200 K and 0.44 eV traps at a temperature of 460 K. 
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Fig. 12. CTI from Analytic Model (AM) including the clock voltage 
effects as a function of temperature for 0.44 eV traps with a concentration of 
10 12 cm -3 and 1% hit (pixel) occupancy at 50 MHz readout frequency in 
comparison with full TCAD simulation results. Two different clock voltages 
(Ym) are shown. 
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Fig. 13. CTI from Analytic Model as a function of clock voltage amplitude 
for both traps 0.17 eV and 0.44 eV with a concentration of 10 12 cm -3 at 
T = 200 K and T = 460 K, respectively, with 1% hit (pixel) occupancy 
and 50 MHz readout frequency. 



VI. Conclusions and Outlook 

Our previous Analytic Models for a CPCCD have been 
extended to include the effect of non-uniform signal shape 
and the effect of realistic clock voltage amplitudes for CTI 
calculations. The signal shape affects the CTI mostly in the 
peak region. A smaller width of the potential well decreases 
the CTI. The inclusion of the clock voltage effects leads 
to smaller CTI values only above the CTI peak position. 
In summary, the Analytic Model has been extended to give 
a more realistic description of the CTI for a CPCCD and 
the results agree better with full TCAD simulations. Overall, 
the Analytic Model predicts well the CTI peak position in 
comparison with a full TCAD simulation. It can produce CTI 
values almost instantly while the full TCAD simulation is very 
CPU intensive. Generally, agreement between the Analytic 
Model and full TCAD simulation results is better for the 0.17 
eV traps than for the 0.44 eV traps. The Analytic Model is 
suited to contribute to future CPCCD developments. 
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